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(1.1) 
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1 .  THE  PROBLEM 

The  equations  studied  for  purposes  of  inversion  were 
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and 


(1.2)  E(r)  =  E‘'(r) 


+  k 


G{r  -  f')  ]  f(f' )  E(f')  df* 


A  rapidly  convergent  algorithm  for  solving  the  forward  problem 
corresponding  to  (1.1)  was  carried  out  in  a. 16  above.  Eg.  (1.1) 
(rsp.  (1.2))  obtains  from  the  Helmholtz  differential  equation 
(rsp.  Maxwell  differential  equations) .  Inversion  involves  the 
reconstruction  of  the  .complex  scalar  function  f(f)  by  applying 
inputs  u< (f )  (rsp.  E^(f))  and  then  measuring  u(f)  (rsp.  E(r)) 
on  the  boundary  of  V. 


2.  REMARK:  In  carrying  out  the  solution  to  these  problems, 
nearly  all  of  our  effort  was  directed  towards  the  inversion  of 
Eg.  (1.1)  for  the  following  reasons: 

(a)  Equation  (1.1)  is  a  scalar  equation,  and  hence  the  amount  of 
work,  the  size  of  the  computer  program  and  storage,  and  the 
amount  of  time  required  to  run  a  program  are  all  considerably 
less  than  the  corresponding  ones  for  Eg.  (1.2). 

(b)  Equation  (1.2)  has  higher  order  singularities  than  Eg.  (1.1), 
and  the  correct  treatment  of  these  requires  additional 
computer  code.  In  particular,  it  requires  the  evaluation  of 
integrals  of  the  type  in  (1.1)  as  well  as  the  evaluation  of 
principal  value-type  integrals.  In  the  references  [a. 12, a. 15] 
above,  we  developed  the  tools  which  are  required  to  handle 
the  principal  value-type  singularities. 

(c)  Except  for  the  problem  size  and  the  evaluation  of  the  princi¬ 
pal  value  integrals,  all  other  aspects  of  the  inversion  of 
(1.2)  are  the  same  as  those  for  (1.1).  That  is,  having 
developed  an  effective  algorithm  for  the  inversion  of  (1.1), 
we  can,  in  effect,  "write  down"  the  corresponding  one  for  the 
inversion  of  (1.2). 


3.  BRIEF  SUMMARY  OF  RESULTS  OBTAINED 

During  the  duration  of  this  proposal,  we  have  developed 
some  effective  algorithms  for  inverting  Eg.  (1.1),  which  are 
described  in  the  papers  [a. 5,  a. 16].  The  other  papers  listed 
above,  which  describe  work  that  was  carried  out  with  at  least 
partial  support  of  this  proposal,  involved  the  development  of  new 
mathematics  and  new  algorithms  that  were  required  for  the 
inversion  of  Eqs.  (1.1)  or  (1.2). 


4.  MORE  DETAILED  SUMMARY  OF  RESULTS  OBTAINED 

We  believe  that  we  have  made  important  advancements  in  developing 
new  methods  for  the  inversion  of  Eqs.  (1.1)  and  (1.2).  We  now 
summarize  some  of  these. 


4.1  Evaluation  and  Inversion  of  the  Laplace  Transform. 

This  is  the  manuscript  a.l.  A  copy  of  this  manuscript  was  sent 
to  your  office  earlier.  We  feel  this  is  a  very  important  paper, 
since  it  makes  the  application  of  Laplace  transform  inversion 
just  as  easy  as  the  inversion  of  Fourier  transforms.  Method  III 
in  this  paper  is  the  simplest,  and  I  expect  that  it  will  become 
a  popular  method  for  Laplace  transform  inversion.  Both  the 
evaluation  and  inversion  of  the  Laplace  transform  are  discussed 
in  a. 10. 

4.2.  Multigrid-Sinc  Methods 

The  combination  of  multigrid  and  Sine  methods  provides  a  natural 
"marriage",  for  solving  both  differential  and  integral  equations, 
as  was  shown  in  a. 3.  We  are  in  the  process  of  developing  an 
effective  computer  algorithm,  which  we  expect  will  make  both 
methods  more  popular. 

4.3.  Sine  Solution  of  Cauchy  Singular  Integral  Equations. 

The  paper  a. 11  consists  of  yet  another  "breakthrough",  since 
since  it  makes  possible  the  very  efficient  (optimal)  solution  of 
a  two-dimensional  problem  via  a  one  dimensional  procedure.  It 
will  thus  make  it  possible  to  develop  methods  for  solving  a  large 
class  of  three  dimensional  elliptic  PDE  problems  via  two 
dimensional  methods. 

4.4.  Poisson  Solver 

We  have  done  long  and  painstaking  work  on  the  paper  b.2,  which  also 
represents  a  "breakthrough"  similar,  yet  different  from  the  one  in 
4.3  above.  This,  too,  will  open  many  new  avenues  for  rapid  solution 
of  partial  differential  equations. 

4.5.  Sine  Package  of  Algorithms. 

We  expect  completion  of  this  work  (b.3  above)  during  the  next 
three  months.  It  will  appear  in  "Sine  Methods  in  Numerical 
Analysis",  a  textbook  which  is  in  the  final  stages  of  completion. 
Many  of  the  procedures  referred  to  above,  as  well  as  others  which 
were  previously  supported  by  ORA  will  appear  in  this  text. 


